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We present a new approximate method for constructing gravitational radiation driven inspirals 
of test-bodies orbiting Kerr black holes. Such orbits can be fully described by a semi-latus rectum p, 
an eccentricity e, and an inclination angle l; or, by an energy E, an angular momentum component 
L z , and a third constant Q. Our scheme uses expressions that are exact (within an adiabatic 
approximation) for the rates of change (p, e, I) as linear combinations of the fluxes (E, L z , Q), 
but uses quadrupole-order formulae for these fluxes. This scheme thus encodes the exact orbital 
dynamics, augmenting it with approximate radiation reaction. Comparing inspiral trajectories, we 
find that this approximation agrees well with numerical results for the special cases of eccentric 
equatorial and circular inclined orbits, far more accurate than corresponding weak-field formulae for 
(p, e, i). We use this technique to study the inspiral of a test-body in inclined, eccentric Kerr orbits. 
Our results should be useful tools for constructing approximate waveforms that can be used to study 
data analysis problems for the future LISA gravitational- wave observatory, in lieu of waveforms from 
more rigorous techniques that are currently under development. 

PACS numbers: 04.30.Db, 04.25.Nx, 95.30.Sf 



I. BACKGROUND AND MOTIVATION 



The capture of stellar-mass compact objects by massive black holes residing in galactic nuclei is expected to be one 
of the most important sources of gravitational radiation for the future LISA space-based detector [QJ|]. Observing 
such events will provide information about stellar dynamics in galactic nuclei, and should make possible precise 
measurements of black holes' masses and spins. Indeed, the waves generated by such a capture will encode a detailed 
description of the black hole's spacetime, making it possible to test whether the "large object" in the galactic nucleus 
is indeed a Kerr black hole as predicted by general relativity, or is some exotic massive compact object |||}. 

Extracting such information will require accurate modeling of the gravitational waveform. The smallness of the 
system's mass ratio (typically, n/M ~ 10~ 4 — 10~ 6 , where /i and M are the masses for the captured body and the 
central hole respectively) allows one to treat the small body as a "test particle" moving in the gravitational field of the 
black hole. In the absence of radiation, the small body moves on a geodesic orbit of the black hole ||. These orbits 
have three integrals of motion (apart from fj,): energy E; angular momentum projected on the hole's spin axis, L z ; 
and Carter's third constant Q, related to the square of the angular momentum projected onto the equatorial plane. 
A body in a generic (eccentric and inclined) Kerr orbit traces an open ellipse precessing about the black hole's spin 
axis, resulting in a complicated overall motion. Astrophysical captured bodies will move in such complicated orbits. 

The integrals of the motion are not constant in the presence of gravitational radiation — they evolve as energy and 
angular momentum are carried away by the waves. Due to the small mass ratio, they should change adiabatically, on 
timescales much longer than any orbital timescale. Hence, the orbit looks geodesic on short timescales. This fact can 
be used to calculate gravitational perturbations induced by the particle at infinity and at the hole's event horizon, 
using the Teukolsky-Sasaki-Nakamura formalism ||. In this way, one can explicitly find the gravitational waveform 
at infinity and compute the corresponding fluxes of E and L z to infinity and into the hole. If one also knows the 
evolution of the Carter constant, then the adiabatic nature of the inspiral allows one to treat the small body's motion 
as an evolution through a sequence of orbits: the body's worldline z(t) is that of a geodesic orbit z geo d(i) whose 
orbital constants are slowly changing: 

z(t) = z g cod[t; E(t), L z (t),Q(t)} . (1.1) 

Computing the inspiral properties is reduced to computing the parameter space trajectory [E{t), L z (t), Q(t)]. 

One can in fact infer the change in Q and thus fix the small body's inspiral in two special cases: orbits that are 
equatorial, and orbits that are circular but inclined. A considerable amount of effort has been devoted to studying 
these orbits and their evolution due to gravitational-wave emission In these special cases, the evolution of 

the Carter constant Q is constrained: it remains constant at Q — (equatorial orbits) or evolves such that the 
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system's eccentricity is constant at e = (circular, inclined orbits) |13|-|l5[. Accurate numerical computations, based 
on extracting E and L z from fluxes of gravitational waves to infinity and down the event horizon, have detailed the 
effects of radiation reaction and the nature of gravitational-wave emission in these cases. 

Unfortunately, this "flux-balancing" prescription fails in general — there is no known method for computing the 
rate Q from the gravitational-wave fluxes in the absence of special constraints. At the moment, the only applicable 
result is a weak- field, quadrupole-order calculation by Ryan |^3| , who used a weak-field radiation reaction force to 
infer Q. Not surprisingly, Ryan's results become increasingly inaccurate and unreliable as the orbit comes closer to 
the black hole. It is likely that a strong-field gravitational self-force prescription will be needed to compute Q. Many 
groups are now working on this problem |ll|. It is generally acknowledged that no result applicable to strong-field 
Kerr orbits should be expected within the next few years. In the meantime, therefore, an investigation of possible 
approximation schemes for describing radiation reaction and wave emission by these orbits is highly desirable. Such 
schemes will play an important role in mapping out the scope of the data analysis task that the LISA community 
faces, making possible a realistic assessment of issues such as the amount of computing power needed, the accuracy 
with which black hole characteristics can be measured, and the difficulty of measuring signals if the inspiral rate is 
large enough to create a confusion-limited background p7| . 

In the remainder of this paper, we present such an approximate scheme. The essential idea is to use the exact Kerr 
black hole geodesies to describe the system's dynamics, but to evolve through a sequence of those geodesies using the 
weak-field quadrupole-order fluxes for E and L z . Because this scheme mixes an exact notion of short-timescale motion 
with an approximate description of the long-timescale radiation effects, we call it a "hybrid" approximation. We find 
that the hybrid approximation faithfully reproduces features seen in the numerical strong-field analyses of radiation 
reaction. For example, we find that the rate of change of eccentricity will typically switch sign prior to plunging into 
the black hole; as a consequence, the orbit has substantial eccentricity near the end of inspiral. Self-consistent leading 
order calculations (which approximate the orbital dynamics as well as the radiation reaction) strongly underestimate 
this residual eccentricity. In some cases, they predict that the orbit is circular at the end of inspiral. This incorrect 
circularization could have a big effect on the waveform models that are used to lay the foundations of LISA data 
analysis, since a circular inspiral produces waves with less interesting harmonic structure than eccentric inspirals. 
We advocate this hybrid scheme as a method that is simple enough to produce waveforms that are "fast and dirty" , 
but accurate enough to qualitatively reproduce features that should exist in real inspirals. We emphasize that the 
hybrid waveforms are not the ultimate models one would want to use as templates for analyzing the LISA datastream. 
Instead, we advocate them as tools for exploring issues in LISA data analysis, as described in the paragraph above. 

The ideas behind the hybrid approximation and key equations are given in Sec. H some of the more cumbersome 



details are presented in Appendix |B|. In Sec. Ill, we then compare this technique's predictions to those of detailed 
numerical calculations for the two special cases that are well-understood now, equatorial orbits and inclined, circular 
orbits. We compare with the leading order results developed by Ryan |l|, and show that the hybrid scheme qualita- 
tively recovers features seen in the strong-field numerical calculations. Our results show that holding the inclination 
angle l constant is more accurate than letting it evolve in the way that the weak-field fluxes "want" it to evolve (as 
compared to strong- field numerical calculations). We argue in Sec. Ill and Appendix [A| that this tells us that the 
"gravitational potential" felt by the inspiraling body is nearly spherical, and argue further that holding l constant 
should work well for arbitrary orbits. 



In Sec. IV, we move to "generic" configurations, studying inspirals through a sequence of inclined, eccentric orbits. 
In most cases, we find that the inspiral trajectories are qualitatively similar to inspiral in the equatorial plane. In 
particular, we find that most configurations plunge into the black hole with substantial residual eccentricity. We also 
map out the range of parameter space where we do not trust the hybrid approximation: when orbits reach too deeply 
into the strong field, or spiral in near inclination 90°, the weak-field fluxes that we use do not appear to be reliable. We 
conjecture in Sec. [v| on how an approximation could be developed to better understand the Carter constant's evolution. 
This approximation is based on the "zoom-whirl" behavior of strong-field eccentric orbits, recently described in Ref. 
|jl2|| . We provide concluding discussion and suggest directions for future work on this problem in Sec. VI. Throughout 
this paper, we use units in which G = c = 1 . 



II. THE HYBRID APPROXIMATION 



Generic Kerr geodesies can be parameterized by a triplet of constant orbital elements: the semi-latus rectum p. the 
eccentricity e, and the inclination angle l. The elements p and e define the orbit's radial turning points, the apastron 
and periastron: 



In the strong field of a Kerr black hole, there are many ways that one could define an "inclination angle" — for 
example, the turning points of the orbit's latitudinal motion, or the angle at which the small body crosses the equator 
as seen by distant observers. We use the following definition: 
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This definition does not correspond to either of these examples, but is very convenient: it depends simply on or- 
bital constants and has a useful intuitive description, suggesting that the Carter constant Q is essentially just the 
square of the angular momentum projected into the equatorial plane. (This description is in fact exactly correct for 
Schwarzschild black holes; for non-zero spin it is not quite correct, but is good enough to be useful. We discuss this 
issue in more detail in Appendix |a[) The orbital elements can be written as functions of (E, L Z ,Q), and vice versa. 
Consequently, we can write their time-derivatives as p = p(p, e, i, E, L Z ,Q), and similarly for e and L. 

As already mentioned, we do not yet know how to accurately calculate Q. It is only known to leading order in M/p 
and in the spin of the black hole |13| . The orbital parameters used in reference [ [13| are an eccentricity e (different 
from e), a semi- major axis a, and an inclination angle i (identical to our i). The two sets of parameters are related by 
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The parameterizations are consistent in the weak field, and are identical for zero spin. Rewriting Ryan's fluxes in 
terms of our parameters yields 
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In the a = limit, Eqs. (2.5) and (|2.6|) reduce to the celebrated Peters-Mathews formulae U} 



The rates g'j = {p, e, i} can be written 



q-j =H 1 {b 3 E + c 3 L z + djQ) 



(2.5) 
(2.6) 
(2.7) 
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The quantities H and bj, Cj, dj are all constructed in a straightforward way from derivatives of E, L Zl Q with respect 
to p, e, the resulting expressions are rather cumbersome, and so are written out in Appendix B. We emphasize that 
these functions encode the exact geodesic motion. 
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The main idea behind the hybrid schem e is simple: calculate the time derivatives <jj using the exact coefficients bj, 
Cj, dj and the approximate fluxes (fe.Ej) - ( [2.7|) . A consistent leading-order calculation (that is, leading order in M/p 
and a/M) would instead approximate the coefficients qj along with the fluxes. Knowing the rates <jj makes it possible 
to build the parameter space trajectories Qj(t) followed by a small body spiraling into a black hole: given initial values 
5j(0), one simply "integrates up" the derivatives dj to generate the inspiral trajectory. For example, a simple-minded 
Euler-method integration would step from parameter space coordinates (i, qj) to (t + St,qj + qjSt). Generalization 
to more sophisticated integration techniques is straightforward. The trajectories qj(t) are the main result of this 
paper. From them, it is a simple matter to compute quantities such as the gravitational waveform generated during 
an inspiral, and thus to begin testing ideas more directly related to data analysis. We will not develop such waveforms 
here, but will instead defer them to a later analysis. 



III. COMPARISON WITH NUMERICAL RESULTS 



The reliability of this method can be assessed by applying it to specific families of orbits where accurate numerical 
results are already known. We first consider equatorial eccentric orbits, recently studied by Glampedakis and Kennefick 
jL2|. Such orbit s always have t = 0° (prograde) ox l — 180° (retrograde), leaving p and e as unspecified parameters. 
Equation (2.14) becomes 
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where H cq — E p L z e — E e L ZtP . The leading-order approximation for these expressions is 
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Note that we could equally well use the corresponding expressions with th e le ading- orde r spin terms included (see 
fl2||), but it turns out that they essentially give the same results as Eqs. ( |3.3|) and (^4). These equations can be 
combined to give a simple expression that describes the orbital evolution on the p — e plane: 



P(e) = Pi 
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where pi and are initial values. We are now ready to compare the inspiral trajectories generated by Eq. (3.5) with 
those obtained by the hybrid scheme. 

Representative results for astrophysically relevant initial parameters are shown in Fig. [j]. We compare the leading- 
order trajectories found using Eq. (3.5) (dotted lines) with the trajectories predicted by the hybrid scheme (solid 
lines). Note that the time dependence of the inspiral is suppressed in this figure: most time is actually spent at large 
p. The total duration of an inspiral scales with M 2 /fj,. The shape of a curve, however, does not depend on this ratio, 
provided that the mass ratio is extreme: these curves are universal for /i <C M. We show inspiral for both prograde 
and retrograde orbits, for black hole spins a = 0.5M and a — 0.9M. 

In all cases, the hybrid and the leading-order calculations agree for p ^> M, as expected. Differences between the 
two methods become apparent in the strong field. The leading-order inspiral trajectory exhibits constantly decreasing 
eccentricity. This is in marked contrast with the rigorous strong-field calculations (numerical and analytical) of Refs. 
0,||JlJ] showing that there exists a region near the separatrix of stable/unstable orbits where e reverses sign: the 
eccentricity should grow near the separatrix. It is very encouraging that the eccentricity does in fact grow when the 
hybrid approximation is used. Moreover, the location of the critical points in these curves where e = is in good 
agreement (at the order of a few percent) with the numerical results of Refs. [p|JT^]; see Table § Three of the four 
cases shown in Fig. |l| appear "good" in the sense that the trajectories appear to agree reasonably well with what we 
expect based on strong-field numerical analyses (cf. Ref. |l^|). The same comparison for the fourth case (a = 0.9M, 
prograde; upper plot in the right hand panel of Fig. Q) reveals that both the eccentricity growth near the separatrix 
and the distance of the critical curve e = from the separatrix are excessive. Progra de o rbits o f ra pidly rotating black 
holes reach rather deep into the black hole's strong field where the weak-field fluxes (|2.5|) and (2.6) cannot be trusted. 



4 



As we shall see when we move on to generic inspirals (Sec. IV), this breakdown of the weak- field flux formulae means 
that the hybrid approximation does not accurately describe the inspiral of shallow inclination orbits (i < 20° or so) 
into rapidly rotating black holes (a > 0.85M or so). 
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FIG. 1. Comparing equatorial inspiral. We show inspiral into a hole with spin a = 0.5Af (left panel) and into a hole 
with a = 0.9A/ (right panel). In each panel, the top half shows prograde inspirals and the bottom retrograde inspirals. In 
each set, the dashed line represents the separatrix separating stable from unstable orbits. We used the hybrid approximation 
discussed in the text to radiatively evolve orbits with initial parameters (pi,ei) — (20A/, 0.4), (20A/, 0.99), (1AM, 0.99), and 
(10Af, 0.99). The inspiral trajectories are shown as the heavy lines in each plot. (The final set is not included in the retrograde 
inspirals since the initial conditions are not stable in those cases.) The dotted trajectories in each plot shows the leading-order 
predictions generated using Eq. (|3.5|). Note the significant qualitative difference between the two calculations at the vicinity of 
each separatrix. Note also the extremely large growth in eccentricity seen in the prograde inspirals for a = 0.9M just before 
reaching the separatrix. Comparison with accurate strong-field numerical results shows that this growth is excessive. 



It is possible to get some insight into the superior qualitative description of the inspiral in the strong field region 
given by our approach. The phenomenon of orbital circularization as a result of some form of dissipation is seen 
in many astrophysical situations, such as that of satellites whose orbits are decaying due to atmospheric friction. 
The reason is that the dissipating mechanism causes the particle to "drop" in its potential well, the usual geometry 
of which ensures that the orbital eccentricity decreases. In our case another mechanism becomes significant as the 
unstable plunge orbit is approached at the end of the inspiral. As this occurs the potential becomes shallower (as 
the minimum turns into a saddle point at plunge), and this tends to increase the eccentricity of the orbit. Shortly 
before plunge this mechanism overcomes the circularizing tendency. It is not surprising that the hybrid approximation 
can qualitatively replicate the eccentricity increasing behavior, because it exactly describes the shape of the orbital 
potential, which is so critical to this effect. 

Table [Il| compares data for p and e. In this sample, the hybrid approach clearly is more accurate than the leading- 
order approximation. This comparison is a very strict test of the accuracy of this scheme. As discussed above, we 
believe that the hybrid approximation is reliable as long as r p > 5M. The weak-field fluxes that we use cannot be 
trusted very deep in the strong field — the spin correction terms in Eqs. (2.5) and ( pii| ) dominate the leading order 
term. The method therefore fails when we push to smaller r p . This effectively constrains the black hole spin to 
a < 0.5M for prograde motion — for larger spins, the innermost stable orbit and hence r p come too close to the 
horizon. For retrograde orbits, the results are much more accurate s ince r p ne ver comes close to the horizon, regardless 
of the spin. Finally, we emphasize the essential role Ryan's fluxes ( |2.5|) and (2.6) play in calculating p and e. Had we 
used instead the Peters-Mathews fluxes, the resulting inspirals would predict a rapid circularization under radiation 
reaction: we find that the Peters-Mathews fluxes reduce the eccentricity to zero well before reaching the saddle point 
of the orbital potential, and so the eccentricity never grows. This is in sharp disagreement with the numerical results. 



5 



A major prediction of the hybrid approximation is that for equatorial orbits the residual eccentricity prior to 
plunge should be substantial, in strong contrast to the prediction of the leading order formula ([3^) . In many cases, 
the leading order results predict that the orbit will actually circularize prior to plunge. Because the harmonic structure 
of a circular inspiral is rather different from that of an inspiral with substantial eccentricity, these results have strong 
implications for the waveform models to be used in LISA's data analysis. 

We next consider circular inclined orbits, which were recently studied by Hughes [fiof . One of the most important 
findings of Ref. JToj] is that the angle i remains almost constant during inspiral, even when the particle is crossing 
strong field regions. 

For these orbits, the rates p and I are given by 

V = H^ C (-L Z , L E + E, L L Z ) , (3.6) 
i = H^ c (L ZtP E - E iP L z ) , (3.7) 

where -ff C irc = E >L L z ^ p — L z ^E tP . In order to obtain these formulae we first expressed Q in terms of E and L z making 
use of the "circular goes to circular" theorems [p^|-p^] ; see Ref. jllj for further discussion. The leading-order expression 
for i is fl 
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p follows from Eq. (|3.3|), setting e = 0. 
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FIG. 2. Comparing circular, inclined inspiral. We show inspiral into a hole with spin a — 0.5M (top) and spin a — 0.9A/ 
(bottom). The solid lines show inspiral using the hybrid approximation; the dotted lines show the leading order inspiral 
prediction. The dashed curve shows the separatrix between stable and unstable orbits. Both approximations show that the 
inclination angle increases, especially right before reaching the separatrix. However, the increase predicted by the leading-order 
prediction is far too large, particularly for rapidly spinning black holes. The inspiral predicted by the hybrid approximation is 
closer to what is seen in rigorously computed inspirals. Nonetheless, it too shows an increase in i that is probably excessive. As 
we argue in the text, holding i constant produces an inspiral sequence that is probably closest of all to strong-field calculations 
and should be acceptably accurate. 
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using the hybrid approximation to the the results obtained using Eqs. (3.3) 
together with accurate numerical results from Ref. [0 . Figure || shows inspirals of circular inclined orbits 
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with our method and using the leading-order formulae. Both approximations predict that i changes in such a way 
as to drive the orbit to an equatorial retrograde configuration (that is, t increases). The two calculations agree at 
large radii. In the strong field, the leading-order formulae break down — the inclination angle tends to increase 
dramatically. The behavior of the hybrid-scheme inspiral is more reasonable. 

Although the hybrid scheme is much better behaved in the strong field, the growth of b we see is still quite a bit 
larger than detailed numerical calculations predict Jl0|Jll[] . Based on those numeric al re sults, a more accurate scheme 
would be to simply require that i remain constant. Combining db/dt = with Eq. ( |2.2| ) yields the rule 

Q=^-L z . (3.9) 

This rule is consistent with our description of Q as roughly the squared component of angular momentum projected 
into the equatorial plane. If the spacetime is perfectly spherical (i.e., Schwarzschild black holes), then Q is exactly 
such an angular momentum component: Qsphcricai = L^, + L? y . It is easy to show that an inspiral in this spacetime 
would proceed at exactly constan t in clination angle: gravitational waves carry off exactly the right amounts of L x 
and L y to hold i constant, so Eq. ( |3.9[ ) is exactly correct in this case. One would expect i to remain nearly constant if 
the spacetime does not deviate too strongly from sphericity. Rigorous numerical results for the circular inclined case 
show that l indeed remains nearly constant; it thus appears that the Kerr metric is not too aspherical over much of 
the inspiral (modulo frame dragging) . Additional evidence for the validity of this statement is given by the discussion 
in Appendix [A]. Since the orbit's eccentricity does not enter this argument at all, it is likely that Eq. (3.9) will work 
well for inclined circular orbits also. 



IV. EVOLVING GENERIC ORBITS 



Having established the reliability and limitations of the hybrid scheme, we move to the main subject of this paper: 
the study of inspirals of test bodies in generic orbits wher e on ly le adin g-order results are currently available p3[ . 

We began this analysis employing Ryan's fluxes, Eqs. (2.5) - fl2.7|) , but quickly faced disappointing results. We 
found that hybrid-scheme inspirals produced with these fluxes did not behave well far from the two limits discussed 
above, particularly in the strong field. For example, the eccentricity tended to grow extremely large very rapidly 
in some cases. The root of the problem lies in the expression for the Q flux, Eq. (2.7), which apparently is not as 
accurate as we would require it to be. The qualitative behavior of our inspirals is more reasonable when the rule given 
by Eq. ( |3.9| ) is used to compute Q instead, forcing i to be constant. Following the discussion at the end of Section III 
and in Appendix [A], it is likely that this rule is accurate enough for our purposes anyway, and so we shall use it from 
this point onward. In all likelihood, detailed self-force calculations will be needed to test the accuracy of the constant 
inclination rule. 



Using the fluxes given in Eqs. ( |2.5| ) and (2.6) with the constant-t rule (3J5) produces inspirals that agree with the 
leading-order results when p ^ M, that smoothly converge to the equatorial case for l — > 0° and i — > 180°, and that 
exhibit an e > region near the separatrix. Examples of this behavior are shown in Fig. || (spin a = 0.5M) and 
Fig. ^ (spin a = 0.9M). Several interesting features can be seen in these figures. The trajectories for b < 90° are 
qualitatively similar to the equatorial, prograde trajectories shown in Fig. [l| In particular, each such trajectory passes 
through a critical point at which e = after which eccentricity grows. The system typically has substantial non-zero 
eccentricity when it reaches the separatrix. Also, note the excessive growth in eccentricity near the separatrix for 
a = 0.9M and t = 30°. At shallow inclination angle, the separatrix is very deep in the black hole's strong field, so the 
inspiral proceeds to small r before plunging. Just as in the case of equatorial orbits for a — 0.9M, the weak-field flux 
formulae that we use cannot be trusted this far into the Kerr black hole strong field. 

The qualitative appearance of the inspirals for l > 90° is quite a bit different from the b < 90° inspirals. Looking 
at the right hand sides of Figs. |^ and [l[ we see that the hybrid approximation predicts that many of these inspirals 
completely circularize prior to reaching the separatrix. We do not believe that this behavior is robust. Indeed, we 
find that the behavior of inspirals exhibits a rather sharp transition as the inclination angle goes from slightly below 
90° to slightly above. This behavior arises from the cost terms in Eqs. (2.5) and (2.6), which switch sign at this 
transition. We thus do not believe that this rapid circularization is physical, but instead attribute it to poor behavior 
of the hybrid approximation at b > 90°. 
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FIG. 3. Comparing generic inspiral at several inclination angles into a hole with spin a — 0.5M. In all plots, the dashed 
line represents the separatrix between stable and unstable orbits; because we force inspiral to lie in surfaces of constant t, 
there exists a separatrix in the p — e plane of each such surface. We use the hybrid inspiral scheme to evolve orbits with 
initial parameters (pi,ei) = (20M, 0.2), (20M,0.4), (20M,0.99), and (14M,0.99), and i = 30°, 60°, 80°, 100°, 120°, and 150°. 
Inspirals that are roughly "prograde" (have t < 90°) are qualitatively similar to the equatorial inspirals shown in Fig. |l|: there 
is a turnaround in the eccentricity evolution, so that all inspirals finish with a substantial non-zero eccentricity. By contrast, 
the roughly "retrograde" inspirals (t > 90°) exhibit rather different behavior: particularly when the inclination is not too far 
from 90°, many inspirals completely circularize, reaching e = 0. As discussed in the text, we believe this behavior is incorrect. 



Having established that the behavior of hybrid approximation inspirals for i > 90° is probably not reliable, it 
is worth re-examining the behavior for i < 90°. Good examples of this beha vior a re t he plots for i = 80° (lower 
leftmost panels of Figs. || and |]). In these cases, the cost terms in Eqs. (2.5) and ( |2.6| ) will be small but positive. 
Indeed, we see that the trajectories shown in this case are somewhat odd. Consider the trajectories that begin at 
(pi,ei) = (20M, 0.2). The points where the eccentricity evolution switches sign are at rather large semi-latus rectum 
(p ~ 16. 5M for a — 0.5M; p ~ 25M for a = 0.9M). This is quite a bit further out than is seen in any analysis of 
radiation reaction on equatorial orbits J7j,||,[l2| . We suspect that this behavior is likewise an artifact of the weak- field 
fluxes, and do not trust the hybrid approximation's predictions for insp iral s near i = 90°. 

W e co nclude that the hybrid inspiral scheme — the weak- field fluxes (2.5) and ( |2.6| ) plus the "constant inclination" 
rule (|3.9| ) applied to exact, strong-field Kerr geodesies — is, in most cases, reliable and accurate enough to be used for 
exploring issues in LISA's data analysis. In some cases, the hybrid scheme will not be accurate enough: the weak-field 
fluxes behave badly when the orbit goes too deep into the strong field, and so we do not trust this scheme's results 
when r p < 5M. Also, the spin correction terms in Eqs. ( |2.5D and (|]^) become either very small or very large relative 
to the leading quadrupole term when l ~ 90°, and so we do not trust the hybrid approximation for inclination angles 
80° < i < 120°. More rigorous strong-field analyses will be needed in order to validate the inspiral behavior at these 
inclination angles. 

In all cases in which the inspiral behavior is reasonable, we find that small body's orbits will have significant 
eccentricity upon reaching the separatrix. Eccentricity will be a significant factor that must be incorporated into 
plans for LISA data analysis. 
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FIG. 4. Comparing generic inspiral at several inclination angles into a hole with spin a = 0.9M. Aside from the faster black 
hole spin, the trajectories shown here have identical initial parameters as those shown in Fig. |j| The inspirals for t < 90° are 
again qualitatively similar to equatorial inspirals, shown in Fig. ^. In particular, we see that at shallow inclination angle, the 
growth in eccentricity is very large. We believe this is because these orbits go so deeply into the strong-field that the weak-field 
flux formulae used in the hybrid approximation are no longer reliable. We also see the rapid circularization of inspirals when l 
is greater than but near 90°, very similar to the behavior encountered for spin a — 0.5M. 

V. CONJECTURE: APPROXIMATE Q FOR GENERIC ZOOM- WHIRL ORBITS 

We have repeatedly emphasized that the hybrid approximation presented here is reliable as long as the orbiting 
body does not come too close to the central black hole. This excludes shallow inclination orbits of rapidly rotating 
black holes — an unfortunate exclusion, since those orbits are in many cases very strongly "stamped" by the features 
of the strong- field spacetime, and thus may be the most interesting inspiral sources observed by LISA j2(J] . These are 
also the orbits for which the "constant inclination Q" rule ( |3.9| ) is most likely to be inaccurate, since they are deepest 
in the Kerr black hole's strong field and are most likely to feel the asphericity of the Kerr spacetime. Ultimately, 
self-force computations will provide the toolsets needed to rigorously compute strong-field inspirals in this regime, 
and will sidestep all difficulties regarding the calculation of Q. In the meantime, while such computations remain 
unavailable, it is worth exploring other possibilities that may provide accurate approximations to the Carter constant's 
evolution. 

Strong-field eccentric equatorial orbits of rapidly rotating holes have a "zoom-whirl" character | fi2[ : orbits near the 
separatrix in the p — e plane "whirl" around the black hole repeatedly near periastron, so that the amount of azimuth 
(j) accumulated in a single radial cycle (apastron to periastron to apastron) is much greater than 2ir. During this whirl 
phase, the body's motion is very close to a circular orbit. Exploratory studies |2lf| show that this behavior carries 
over to non-equatorial orbits, particularly for fairly shallow inclination angle. 

The equatorial zoom- whirl orbits studied in Ref. [ fl2[ radiate energy and angular momentum as if they were nearly 
circular: 

E « ^L z , (5.1) 

where fi^ is the frequency associated with azimuthal motion. This property follows quite naturally from the motion 
of a test-body in a zoom- whirl orbit: a large fraction of the orbital period is spent "whirling" in the vicinity of the 
periastron, where the motion is nearly circular. This is also the part of the orbit where the body is closest to the black 
hole and thus where most of the radiation will be generated. Thus, the radiation from a zoom-whirl orbit should be 
very similar to radiation from a circular orbit, which is exactly what Eq. (^T]) suggests. Extrapolating this behavior 
to generic zoom- whirl orbits, we expect that most of the radiated E, L z , and Q will come from the motion of the 
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body near a (generalized) separatrix in the (p, i, e) phase space. The whirl motion of such orbits will be well-described 
as nearly circular and inclined. 

Following Kennefick and Ori jL4| , we can write the Carter constant as 

Q = G{r,E,L z )-Au 2 r7 (5.2) 

where 

G = A" 1 [£(r 2 + a 2 ) - aL z } 2 - (L z - aE) 2 - r 2 , (5.3) 

with A = r 2 — 2Mr + a 2 and u r denoting the radial component of the body's four- velocity. (The function that we 
have denoted G is written H in [fu}.) It is then straightforward to show that fl4|| 

2Y,y r 

Q = G.eE + G, l _L z —F r , (5.4) 

u l 

where £ = r 2 + a 2 cos 2 8 and F r is the radial component of the self force. It is the unknown last term in this equation 
that presently prohibits the calculation of Q for generic orbits. For strictly circular orbits, on the other hand, this 
term is absent since u r = 0. The remaining expression Q = Q(E, L z ) allows the explicit calculation of Q; this is how 
Hughes evolves circular, inclined orbits by reading the fluxes E and L z at infinity and down the hole ]Tb] , pj| . 

For a zoom- whirl orbit and for motion near the periastron, r w r p , so we should have u r w 0; consequently, 



the unknown term in (5.4) should be negligible. Our conjecture is that the resulting expression for Q describes the 
evolution of the Carter constant for all generic zoom-whirl orbits and with increasing accuracy as the orbit approaches 
the separatrix. We emphasize that this approximation should hold even for orbits deep in the black hole's strong- 
field. This conjecture could become a practical tool once a cod e th at calculates E and L z for generic orbits is 



developed. Furthermore, a direct comparison between (iS.) and (5^) should be a useful guide for the accuracy of 



the l = constant rule in strong-field situations. Future computation of the self force will provide the ultimate test for 
both approximations. 



VI. CONCLUDING DISCUSSION 



The hybrid approximation presented in this paper should be a valuable tool for the gravitational-wave astrophysics 
community as studies of LISA's data analysis requirements begin, and thence models of the waves generated by com- 
pact bodies spiraling into massive black holes become needed. Such approximate "kludged" waveforms are obviously 
too crude to actually be used in future measurements of compact bodies spiraling into massive black holes; data 
analysis strategies based on waveforms built from rigorous strong-field radiation reaction will be needed. Waveforms 
from approximate inspiral models should be adequate to begin the process of developing a data analysis infrastructure. 
For example, they will be useful for counting the number of analysis filters needed, assessing the computational cost 
of data analysis, and experimenting with data analysis techniques. As rigorous and reliable waveform models become 
available, they can simply be dropped into the codes and infrastructure developed using the hybrid approximation. 

Because this approximation combines the exact strong-field Kerr geodesies with weak-field radiation reaction for- 
mulae, it is somewhat limited: inspiral cannot go too deeply into the strong field, thereby making it inaccurate for 
shallow (l < 20°) inspirals of rapidly rotating (a > 0.85A/ ) holes. Also, the cos t dependence of terms within the flux 
formulae behaves badly near t ~ 90°, so that the approximation is probably not reliable within an inclination range 
80° < i < 120°. Despite these limitations, we have found the hybrid approximation reliably and robustly duplicates 
many of the inspiral properties that we expect will be found when it is possible to study these systems using truly 
strong-field gravitational radiation reaction. In particular, it produces inspiral trajectories that retain substantial 
non-zero eccentricity just before plunging into the hole, as is expected from strong-field analyses in special cases 
J^,[s|JT2]] . We emphasize this point because the harmonic structure of gravitational waves from eccentric orbits is quite 
a bit different from that of waves generated by circular orbits. The residual eccentricity of typical inspirals is likely 
to impact data analysis rather strongly. 

Obviously, waveforms constructed from hybrid approximation inspirals are by no means the ultimate models that will 
be needed for LISA data analysis — we strongly advocate continuing to develop techniques for understanding strong- 
field radiation reaction. Future insight from such studies may make it possible to improve the hybrid approximation. 
Even when strong-field radiation reaction is mature enough to model arbitrary compact body inspirals, the calculation 
may be computationally expensive, so that an approximation scheme of some sort may remain useful. 

Although our overall goal is to provide a tool that can be used to model the gravitational waves produced by 
compact body inspiral, we have presented no such waves in this paper. That is the next step. The calculations 
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we have presented explicitly construct the parameter space trajectories [E(t),L z (t),Q(t)] describing an inspiral. It 
is then a simple matter to combine such a trajectory with the geodesic equations for the Kerr spacetime f22]| to 
produce the Boyer-Lindquist coordinate space inspiral [r(t), 9(t), 4>{t)]. This set of functions specifies the worldline of 
the inspiraling body, and one can use it to compute the gravitational waveform seen by distant observers (see, for 
example, Ref. ^3|). Codes to perform this next step are under development results should be presented in the 
near future. 
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APPENDIX A: 



The occurrence of a "third" orbital constant Q in axisymmetric gravitational fields is not an exclusive feature of 
general relativity. For example, it is familiar from Newtonian celestial mechanics applied to orbital motion in galactic 



gravitational potentials (see, for example Ref. |24|, where the third constant is denoted /). The departure of Q from 

-2 

J y 



L 2 + L 2 can then be attributed to the asphericity of the potential. If such a potential does not deviate very much 



from sphericity, L 2 (the square of the total angular momentum) turns out to be almost constant, so that Q should 
be, after all, nearly L 2 — L 2 . 

It is straightforward to check whether this behavior of L 2 occurs in Kerr spacetime. The definition we use for L 2 
is identical to that used in Schwarzschild spacetime, 

L 2 =p 2 + ( S m9)- 2 p 2 , (Al) 

where p a denotes the particle's four-momentum. For the Carter constant Q we have |22| 

Q=p 2 + (sm9y 2 p 2 -p 2 +a 2 co S 2 d(l-E 2 ) . (A2) 

Combining these two expressions gives 

Q = L 2 - L 2 Z + a 2 cos 2 0(1 - E 2 ) . (A3) 

In other words, Q can be interpreted as the projection of the total angular momentum on the equatorial plane, modulo 
the "aspherical" term a 2 cos 2 9(1 — E 2 ). This interpretation makes sense when the aspherical term is small — that is, 
when a <C M (slow rotation) and/or E « 1 (weak-field orbits). In practice, we find that this term is often significantly 
smaller than the preceding terms even for motion in strong-field regions of rapidly rotating holes. We illustrate this 
in Fig. |^, showing how the quantity SL 2 = L 2 / (Q + L 2 ) — 1 varies with time for a variety of generic orbits around a 
rapidly spinning hole. 

Examining Fig. ||, we see that L 2 deviates very little from Q + L 2 even when the small body is deep in the black 
hole's strong field — in this sample, the difference is no more than about 1%. This shows that interpreting Q as a 
squared projection of angular momentum into the equatorial plane is sensible. Because Q + L 2 is a constant quantity, 
this figure also demonstrates that L 2 is nearly constant. This is exactly what we expect for motion in an axisymmetric 
potential that is almost spherical. These pieces of evidence suggest that the Kerr spacetime is not as "aspherical" 
as we might have expected, at least for the purposes of this argument, lending credence to our suggestion that the 
"i = constant" assumption should be reliable, as discussed in the paper's main body. 
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FIG. 5. Examining our notion of "total angular momentum" for strong-field Kerr black hole orbits. Each panel compares 
the angular momentum squared L 2 = pg + (sin#) -2 p^, to Q + L 2 : the quantity plotted is 8L 2 = L 2 /(Q + L 2 ) — 1. The top 
panel shows these quantities over an orbit with p — 10M, the center panel the quantities over an orbit with p = 5M, and the 
bottom over an orbit p — AM. In all cases, the orbits have eccentricity e = 0.7, inclination i = 60°, and are about a hole with 
spin a — 0.999M. Even deep in the strong field, L 2 differs very little from Q + L 2 — the greatest deviation in this sample is 
about 1%. Since Q + L 2 is a constant by definition, this also shows that L 2 is approximately conserved over the orbit. 



APPENDIX B: 



This appendix contains explic it expressions (in terms of E, L z , Q and their derivatives) for the various functions 
appearing in the formulae (2.14) for the rates p, e, I. First, 



H — Q ,pE^ e L z ^ L — Q tP E^L z ^ e — Q , e E iP L Zt , + Q . e E^,L z . p + Q.i.E. p L z . e — Q ti E <e L ZtP 



(Bl) 



For p we find 



For e, 



bp 


= Q,uL z _ e 


— Q,eL z 




= E^ L Q : e - 




dp 


= E^ e L z ^ 


— E,lL z , 


be 


L Zj t,Q^p 


~ Q,iL Zt 




= Q,lE <p - 


- E, L Q^p 


d e 


= E^L Z}P 


— E tP L z , 



(B2) 
(B3) 
(B4) 



(B5) 
(B6) 
(B7) 



Finally, for I the coefficients are, 
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Q ,p^z,e 




— Q,p^,e 


~ Q,e^.p i 


d, 


= E^ p L z ^ e 





(B8) 
(B9) 
(BIO) 



a/M 


e 


Pcrit/M (numerical) 


Pait/M (hybrid approx.) 


Fractional difference 





0.2 


6.76 


6.92 


0.0237 





0.4 


6.99 


7.13 


0.0200 


0.5 


0.3 


4.85 


5.06 


0.0433 


0.5 


0.5 


5.08 


5.21 


0.0250 


-0.99 


0.3 


10.25 


10.53 


0.0273 


-0.99 


0.5 


10.59 


10.78 


0.0179 



TABLE I. Comparing critical curve values p cr it for equatorial eccentric orbits. These are the p values at 
which the eccentricity evolution switches sign, beginning to grow rather than shrink. We show the values 
of Pcrit calculated numerically [||,[t2| (third column) and using the hybrid approximation (fourth column), 
for a variety of black hole spins (negative a/M represents retrograde orbits) and eccentricities. The fifth 
column shows the fractional difference between the numerical and the approximate results, (approximate 
-numerical) / (numerical) . 



Q/iVl 


^ / A/7" 
p/M 


e 


Calculation 


[lVl/fl)p 




Frac. diff. in p 






Frac. diff. in e 





7.505 


0.189 


Numerical 


-7.475 x 10" 


-12 




-1.967 x 10" 


-a 










Hybrid 


-6.859 x 10" 


-2 


0.0824 


-1.291 x 10" 


■3 


0.3434 








Leading order 


-2.957 x 10" 


■2 


0.6044 


-1.159 x 10" 


■3 


0.4108 





6.9 


0.4 


Numerical 


-4.240 x 10" 


-I 




+1.047 x 10" 


-a 










Hybrid 


-3.056 x 10" 


-1 


0.2792 


+1.506 x 10" 


-2 


-0.4384 








Leading order 


-3.420 x 10" 


-2 


0.9193 


-2.929 x 10" 


3 


1.2797 


0.5 


6.5 


0.4 


Numerical 


-5.999 x 10" 


-> 




-5.155 x 10" 


•i 










Hybrid 


-4.606 x 10" 


-2 


0.2322 


-3.356 x 10" 


■3 


0.3490 








Leading order 


-4.091 x 10" 


-2 


0.3181 


-3.719 x 10" 


3 


0.2786 


0.5 


15 


0.4 


Numerical 


-3.371 x 10" 


3 




-1.341 x 10" 


-4 










Hybrid 


-3.358 x 10" 


-3 


0.0039 


-1.334 x 10" 


-4 


0.0052 








Leading order 


-3.328 x 10" 


3 


0.0128 


-1.311 x 10" 


-4 


0.0224 


0.5 


4.8 


0.3 


Numerical 


-6.354 x 10" 


-1 




+1.369 x 10" 


"1 










Hybrid 


-4.858 x 10" 


-1 


0.2354 


+3.519 x 10" 


-2 


-1.5705 








Leading order 


-4.849 x 10" 


■2 


0.9237 


-4.432 x 10" 


■3 


1.3237 


0.9 


5 


0.4 


Numerical 


-7.507 x 10" 


-> 




-9.266 x 10" 


■a 










Hybrid 


-4.617 x 10" 


-2 


0.3850 


-1.965 x 10" 


■3 


0.7879 








Leading order 


-2.698 x 10" 


-3 


0.9641 


-9.732 x 10" 


-5 


0.9895 


-0.5 


10 


0.4 


Numerical 


-3.115 x 10" 


-12 




-1.379 x 10" 


-a 










Hybrid 


-2.494 x 10" 


-2 


0.1993 


-9.107 x 10" 


-4 


0.3396 








Leading order 


-1.337 x 10" 


-2 


0.5708 


-7.931 x 10" 


-4 


0.4249 


-0.99 


10.5 


0.4 


Numerical 


-7.516 x 10" 


-> 




-5.312 x 10" 


-4 










Hybrid 


-5.506 x 10" 


-2 


0.2674 


-5.223 x 10" 


-4 


0.0168 








Leading order 


-9.704 x 10" 


-3 


0.8709 


-5.461 x 10" 


-4 


0.0281 


-0.99 


15 


0.4 


Numerical 


-5.766 x 10" 


3 




-2.141 x 10" 


-4 










Hybrid 


-5.295 x 10" 


-3 


0.0817 


-1.875 x 10" 


-4 


0.1242 








Leading order 


-3.328 x 10" 


-3 


0.4228 


-1.311 x 10" 


-4 


0.3877 



TABLE II. Comparing the rates p, e for several equatorial eccentric orbits. The fifth column in this table 
shows p; the seventh column shows e. Within each section of the table, the first row of columns five and seven 
contains accurate numerical data from |8|Jl2[, the second row shows data using the hybrid scheme outlined 
in this paper, and the third row shows data using quadrupole order results. The sixth and eighth columns 
show the fractional differences between the two approximation schemes and the accurate numerical results. 
In all cases but one, the hybrid approximation is closer to the accurate numerical calculation, sometimes 
substantially so. 
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a/M 


p/M 


i(degrees) 


Calculation 


(M/n)p 




Frac. diff. in p 


(M 2 /li)L 




Frac. diff. in I 


0.95 


7 


62.43 


Numerical 


-4.657 x 10" 


-2 


— 


1.207 x 10" 


-4 


— 








Hybrid 


-4.497 x 10" 


- 2 


0.0344 


2.639 x 10" 


-4 


1.1864 








Leading order 


-2.750 x 10" 




0.4095 


3.080 x 10" 




1.5518 


0.05 


7 


60.17 


Numerical 


-1.096 x 10" 


- 1 


— 


1.087 x 10" 


- b 


— 








Hybrid 


-1.048 x 10" 


- 1 


0.0438 


1.207 x 10" 


-5 


0.1104 








Leading order 


-3.676 x 10" 


. 2 


0.6642 


1.587 x 10" 




0.4500 


0.5 


10 


67.56 


Numerical 


-1.583 x 10" 


-2 


— 


1.546 x 10" 


- b 










Hybrid 


-1.645 x 10" 


-2 


0.0392 


2.043 x 10" 


-5 


0.3215 








Leading order 


-1.194 x 10" 




0.2457 


2.377 x 10" 




0.5375 


0.5 


10 


126.76 


Numerical 


-2.329 x 10" 


-2 


— 


1.892 x 10" 


-b 


— 








Hybrid 


-2.341 x 10" 


-2 


0.0051 


1.643 x 10" 


-5 


0.1316 








Leading order 


-1.414 x 10" 




0.3929 


2.060 x 10" 




0.0888 


0.9 


10 


74.07 


Numerical 


-1.544 x 10" 


-2 


— 


2.715 x 10" 


- b 


— 








Hybrid 


-1.567 x 10" 


-2 


0.0149 


3.857 x 10" 


-5 


0.4206 








Leading order 


-1.169 x 10" 




0.2429 


4.452 x 10" 




0.6398 


0.9 


10 


131.57 


Numerical 


-3.253 x 10" 


-2 


— 


3.887 x 10" 


- b 


— 








Hybrid 


-3.082 x 10" 


-2 


0.0526 


2.612 x 10" 


-5 


0.3280 








Leading order 


-1.545 x 10" 




0.5250 


3.464 x 10" 




0.1088 


0.5 


6 


48.33 


Numerical 


-1.237 x 10" 


- 1 


— 


1.410 x 10" 


-4 


— 








Hybrid 


-1.135 x 10" 


- 1 


0.0824 


2.614 x 10" 


-4 


0.8539 








Leading order 


-4.440 x 10" 




0.6411 


3.190 x 10" 




1.2624 


0.5 


6 


67.81 


Numerical 


-2.020 x 10" 


- 1 


— 


2.094 x 10" 


-4 


— 








Hybrid 


-1.779 x 10" 


- 1 


0.1193 


2.992 x 10" 


-4 


0.4288 








Leading order 


-5.082 x 10" 


2 


0.7484 


3.954 x 10" 


-4 


0.8882 


0.9 


6 


54.64 


Numerical 


-7.846 x 10" 


-2 




2.015 x 10" 


-4 










Hybrid 


-6.950 x 10" 


2 


0.1142 


5.486 x 10" 


-4 


1.7226 








Leading order 


-3.598 x 10" 


2 


0.5674 


6.268 x 10" 


-4 


2.1107 


0.9 


6 


99.55 


Numerical 


-74.32 






6.337 x 10" 


-4 










Hybrid 


-48.02 




0.3539 


5.241 x 10" 


-4 


0.1729 








Leading order 


-6.593 x 10" 


2 


0.9991 


7.580 x 10" 


-4 


0.1961 



TABLE III. Comparing the rates p, I for several inclined circular orbits. The fifth column shows p; the 
seventh column shows L. Within each section of the table, the first row of columns five and seven contains 
accurate numerical data from jjl]] , the second row shows data using the hybrid scheme outlined in this paper, 
and the third row shows data using leading order results. The sixth and eighth columns show the fractional 
differences between the two approximation schemes and the accurate numerical results. In most cases in this 
sample, the hybrid scheme performs much better than the leading-order approximation when compared to 
the rigorous numerical data. The only case which this is not true is for L of retrograde orbits. Nevertheless, 
this small inaccuracy has no impact on the calculation of generic inspirals, as we assume that i — constant. 
Note the enormous difference between the numerical and the leading order results in the Table's final entry. 
This is because that point is fairly close to the separatrix between stable and unstable orbits. Since the 
leading-order calculation has no notion of this separatrix, it is particularly inaccurate here. 
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